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Neutron stars possess some of the strongest magnetic fields known in the universe. The surface magnetic 
fields of radio pulsars are estimated to be in the range 10** to 10" Gauss, with 10'" Gauss being the typical 
value. Magnetars, a class of neutron stars with even stronger magnetic fields, ~ 10'^ Gauss, are believed to be 
"magnetically powered" stars, deriving most of their radiative luminosity at the cost of their magnetic fields. 
The origin of the strong magnetic fields of neutron stars, in particular those of magnetars, has essentially been 
an open question for decades. In this paper we explore the possibility that a magnetar may owe its strong field 
to a magnetized core which, as indicated by certain equations of state, may form due to phase transitions at high 
density mediated by strong interaction within a sufficiently massive neutron star. We argue that the field derived 
from such a core could explain several inferred evolutionary behaviors of magnetars. 
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I. INTRODUCTION 

Magnetars are neutron stars with surface magnetic fields a 
thousand times larger than that of an average pulsar. However, 
this is not their most distinctive property. While in pulsars, the 
loss of rotational energy can account for all of the observed 
radiation, in magnetars this is not so. The magnetars exhibit 
spin down ages ~ 10"* years. Over this period, they emit a qui- 
escent radiative luminosity ~ 10"'^ - lO-'^ erg/s. In addition, 
some of them emit repeated flares or bursts of energy typi- 
cally ~ lO'^^-lO"*'* erg, and at times much higher: ~ 10'"' erg 
(e.g. the very energetic flare of SGRl 806-20 on 27 December 
2004 illll). The energy emitted in both quiescent emission 
and flares far exceeds the loss in their rotational energy over 
the same period. The only known source of energy for these 
emissions is their magnetic energy IHH], yet their spin his- 
tory appears to present no clear evidence of a decrease in their 
surface magnetic fields with time lUt] . 

The range of spin periods in which magnetars have so far 
been found is surprisingly narrow, between 5 and 12 seconds. 
While the absence of shorter periods can be attributed to rather 
fast spin down due to their strong magnetic field, the upper 
cutoff in the spin period distribution implies a sudden ceasing 
of magnetar activity while their magnetic fields are still strong. 

The distinct radiative property of magnetars indicates that 
they have an internal state different from that of ordinary ra- 
dio pulsars. During their active period magnetars must have 
a transient store of free energy which dissipates as the ob- 
served radiation. In this paper we show that if in the collapse 
of the progenitor to a neutron star, the strong interaction cre- 
ates a strongly magnetic core, most of the magnetar behavior 
can be explained. The magnetic field of the core will initially 
be screened by the conducting exterior. As this configuration 
relaxes, magnetic free energy is released and drives the mag- 
netar activity. In section below, we discuss the evolution 
of the magnetic field expected in the presence of a strongly 
magnetised core. In section |III] we discuss how, in certain 
equations of state of neutron star matter, such a strongly mag- 
netized core could naturally arise. Section|IV]summarises our 
conclusions. 



II. MAGNETIC FIELD EVOLUTION IN MAGNETARS 

We start with the assumption that in the collapse of the pro- 
tostar, a magnetized core is created in the newly formed neu- 
tron star with typical fields at the core surface of ~ 10'^ - 
10'^ gauss. Any dynamically generated field due to spin align- 
ment via strong interaction is expected to be of this magnitude, 
which roughly equals the product of the nuclear magneton and 
the baryon number density in the neutron star interior 

It is unlikely that the magnetic moment of the core will 
manifest itself in a corresponding magnetic field on the neu- 
tron star surface immediately upon its formation. The mate- 
rial surrounding the core is a pre-existing, highly conducting, 
neutron-proton-electron (npe) plasma. As the core undergoes 
a phase transition to the magnetic phase, this suiTounding npe 
plasma will see the magnetic field rise at the core boundary 
in a time scale set by causality, Jrcore/c ~ 10"^ sec, which 
is much longer than its inverse plasma frequency, ~ 10"^' s. 
Here rcore is the radius of the magnetic core. In reality the rise 
time of the field at the core boundary may be much longer: 
the phase transition may start in small islands, which will then 
have to come together and merge to form a coherent macro- 
scopic structure. The islands may in principle have arbitrary 
magnetic orientation with respect to each other, but any pre- 
existing magnetic field, inherited by the star from its progen- 
itor, will serve to define the orientation of the magnetic mo- 
ments of these small regions, and hence also determine the 
orientation of the net magnetic moment of the core. It there- 
fore appears that the core magnetic field will establish itself 
over a timescale long enough for the suiTounding npe plasma 
to set up screening currents, as per Lenz's law, to prevent the 
penetration of the magnetic flux through it. In this state, the 
magnetic field due to the core will be squeezed into a tight 
volume inside the star, around rcme, while an outside observer 
will have no knowledge of it. This magnetic configuration has 
a much higher energy than a relaxed, dipole distribution of the 
magnetic field outside the core. Typically, if the magnetic flux 
is confined within a radius r, its magnetic energy would be 
proportional to r"'. For a core radius $50% of the stellar ra- 
dius, the excess magnetic energy would be of the same order 
as the energy in the relaxed dipole configuration. 
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With time, the magnetic field will tend to relax to the dipole 
configuration of lower energy, releasing the excess magnetic 
energy in the process. This will happen because the screening 
currents have a finite lifetime-they will be eventually killed by 
dissipativenrocesses such as ohmic dissipation and ambipolar 
diff'usion Oa]. Time scales for these processes, relevant to npe 
plasma in the interior of a neutron star, have been worked out 
by Goldreich and Reisenegger [7]. Their estimates show that 
for typical temperatures in the neutron star interior, the ohmic 
dissipation time scale would be very long: ilO" y. Ambipo- 
lar diffusion will, however, play a very important role, with a 
dissipation time scale ~ 5 x 10'^ / {B^^T^ ^) y, where Z?i6 is the 
local magnetic field strength in unit of 10'^ G and Tg 5 is the 
temperature in units of 10^^ K, a typical value in the interior 
of a young neutron star This would suggest that the strong 
magnetic field will emerge at the neutron star surface in the 
time scale of a few times lO** y. The excess magnetic energy 
will be dissipated over this period, driving an average radia- 
tive luminosity ~ lO-'^ - lO-'^ erg/s. These estimates are in 
good agreement with the inferred ages and observed average 
luminosities of magnetars. For a radius of the core about half 
the stellar radius or less, the total free energy stored in the ini- 
tial magnetic configuration would be sufficient to explain the 
observed level of magnetar activity over a duration of J 10"^ yr. 
The ongoing dissipation will also serve to keep the stellar inte- 
rior hot enough for ambipolar diffusion to proceed in a reason- 
ably short time scale, despite its extreme temperature sensitiv- 
ity. If the local temperature tends to rise significantly above 
10**^ K, cooling via neutrinos will be dominant and will main- 
tain the temperature at this level. This may well happen in the 
early stages, when the energy is released deep into the stellar 
interior |8]. In this phase the X-ray emission from the surface 
may be relatively faint. The young (~ 1700 y), high magnetic 
field (~ 4 X 10^^ G) radio pulsai- 111 19-6127 with faint X- 
ray emission from an unusually hot surface may in fact be an 
example of a star passing through this phase |9]. 

As the excess magnetic energy is dissipated, the magnetic 
field at the stellar surface will rise, eventually reaching the 
full dipole value corresponding to the core magnetic moment. 
This is contrary to what has been normally assumed in the lit- 
erature - that the magnetic field of a magnetar should decrease 
with time as a consequence of dissipation of magnetic energy. 
One way of estimating the magnetic fields of magnetars is 
from their spin-down rate, attributing the spindown torque to 
that due to vacuum magnetic dipole radiation or equivalent. 
The same technique is used to estimate the magnetic fields of 
radio pulsars |10i]. The recorded spin-down history of mag- 
netars, however, shows no evidence of a secular decrease in 
their surface magnetic field. If at all, there is a mild evidence 
in favor of the magnetic field strength increasing with time [5J. 
A similar trend of increasing field strength has been noted in 
certain glitching radio pulsars too ifTlll . 

Magnetars can be classified into two broad categories: 
Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars 
(AXPs) (see Woods and Thompson II12I1 for a recent review). 
SGRs are objects that show relatively frequent burst and flare 
activity. They also possess, on an average, somewhat stronger 
magnetic fields. AXPs represent a quieter magnetar popula- 
tion. Burst activity in them is rare, and their average field 
strength is smaller than that of SGRs. In the conventional 



scenario, AXPs would be the later evolutionary products of 
SGRs 1 12]. In our model, SGRs, with stronger fields, would 
be older than AXPs. As these neutron stars are born in super- 
nova explosions, and happen to be just ~ 10"* y old, they are 
expected to be found associated with the remnants of these 
supernovae. Association with Supernova Remnants (SNRs) 
has indeed been found in several of these cases. Interestingly, 
the AXPs are found to be located close to the centers of the 
remnants while SGRs are often located far from the centers 
of the associated SNRs |5]. Given that all neutron stars re- 
ceive some kick at birth, these observed offsets with respect 
to the SNR centers could be interpreted as the SGRs being 
older, and hence having had a longer time to drift away |5]. 
SGRs have smaller spindown ages (= P/2P) than AXPs, but 
as noted by Thompson, Lyutikov and Kulkarni |5] this may 
simply reflect an accelerated growth of the surface magnetic 
field in this phase, which is expected to happen as the strong 
magnetic field breaks through the upper crust |3, 4]. 

In the time sequence leading to the emergence of the mag- 
netic field at the surface of the star, the leading ambipolar dif- 
fusion zone moves out from the interior of the star to the outer 
surface. Younger objects, where the activity is deeper sub- 
surface, are expected to have a relatively quiet existence, the 
dissipated heat being conducted out to the surface and radi- 
ated away, powering a quiescent X-ray luminosity. As the ob- 
ject ages and the strong field moves through the outer crust, 
mechanical disturbances of the crust could be triggered by 
the magnetic pressure, leading to glitches and flares. This 
trend would qualitatively fit the difference in the level of burst 
activity observed in AXPs and SGRs. The yield strength 
of the upper crust would be unable to support stresses for 
magnetic field gradient across the crust of ^lO'"' G, as the 
maximum stress that the crust can support is estimated to be 
SIO^^ dyne cm"^ 1 13]. It is therefore likely that in this phase of 
evolution the magnetic field will eventually attain its relaxed 
configuration via several hundreds of SGR-type bursts. 

We have mentioned above that the upper cutoff of the spin 
period distribution of magnetars at ~ 12 s implies a sudden 
ceasing of magnetar activity while their magnetic fields are 
still strong. For this there has so far been no natural expla- 
nation (see 1 14] for a discussion). In our model, this abrupt 
stoppage of magnetar activity is in fact expected. Once all 
the screening currents dissipate away and the magnetic field 
reaches its fully relaxed configuration, there is no more free 
energy available in the magnetic field to power the magnetar 
activity. Upon reaching this state, the luminosity of the mag- 
netar will plummet and it will disappear from view. As the 
time scale for the emergence of the field is inversely propor- 
tional to and the spindown torque on the star is proportional 
to B^, it would be natural to expect that at the end of the mag- 
netar activity the spin periods of these stars will be roughly 
similar 1 15]. 

In the post-magnetar phase the spin periods of these neu- 
tron stars will continue to lengthen, and they may function as 
long period rotation-powered radio pulsars. Careful surveys 
sensitive to pulsars with tens of seconds period may be able to 
find them. It is to be noted, however, that the active lifetimes 
of pulsars is inversely proportional to their magnetic field and 
hence the population of such long period, high field pulsars 
is expected to be much smaller than the garden variety radio 
pulsars with teragauss surface field. 
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III. ORIGIN OF MAGNETIC CORES IN NEUTRON STARS 

In the discussion above we have assumed the presence of a 
strongly magnetized core in the neutron stars that act as mag- 
netars. There are several scenarios in which such core mag- 
netic fields can be created dynamically via a phase transition 
at high density for strongly interacting matter To name a cou- 
ple of possibilities, 

1. In classic old fashioned nuclear high density ground 
states, large macroscopic magnetic moment can be gen- 
erated by the alignment of nuclear spins, often accom- 
panied by pion condensation (e.g. ImflTll ). 

2. It is considered more plausible that at typical neutron 
star central densities (~ 5 - 6 times nuclear density) we 
expect the nucleons to dissolve into quark matter We 
have found that the likely ground state of such matter 
at these densities is a neutral pion condensate lfl8[[l9ll . 
Such a state carries a large magnetic moment and hence 
can give rise to a strongly magnetized core. 

Both these cases involve a phase transition to a condensate, 
with condensation energies of order several tens of MeV per 
baryon (see, e.g., ref. [201), corresponding to transition tem- 
peratures in excess of lO'^ K. The interior temperature of a 
neutron star would fall well below this within a day after for- 
mation ll2lll22ll . resulting in the formation of a magnetic core. 

Several other equations of state of quark matter have been 
considered in the literature. One that has attracted much in- 
terest of late is that with a diquark superconductivity [23, 23). 
A transition to such a state is expected at high density (~ 10 
times nuclear density). Ferrer and de la Incera ll25ll have con- 
sidered a scenario where strong magnetic fields can be pro- 
duced by the generation of gluon vortices in the colour super- 
conducing cores of neutron stars. This mechanism achieves 
an amplification of a pre-existing strong (~ 10'^ G) field. 
On the other hand, the magnetic pion condensate mentioned 
above can spontaneously generate magnetar strength fields 
(see fig. 1). 

We now present some considerations that will demarcate 
regions in the parameter space that can support quark matter 
Pion Condensed cores and those which are purely nuclear 

The magnetic pion condensed ground state of quark matter 
from an effective intermediate chiral lagrangian was first con- 
sidered by Kutschera et al |19] for two-flavor quark matter. 
We have generalised the effective lagrangian, L, to 3 flavours, 
which can be used over a large range of densities, from nuclear 
matter to quark matter. We have also, for the first time, fixed 
all the coupling parameters of, L, from low energy hadronic 
data (see ll26ll for details). One of the parameters is the (tree- 
level) mass of the Sigma particle which is related to nn scat- 
tering, and is estimated to be around nia- ~ 800 MeV |27]. 

We have carried out this calculation for charge neutral, 
three-flavor quark matter in /3 equilibrium. This includes 
strange quarks as well as one-gluon exchange interactions. 
We have shown in an earlier work jlSl] that if mg. exceeds 
~ 700 MeV, then according to this mean field theory, the 
quark matter, rather than being in the conventional strange 
quark matter (SQM) state |28], would find it energetically 
favorable to organize itself into a chiral n'^ Pion Condensate 
(PC) which is spin polarized. Our calculations ll26ll show that 



in the PC state that is realized, at the densities of interest, in 
the presence of the 7r° condensate the strong interaction causes 
the u-quark quasiparticles to align their spins along the con- 
densate wave vector q and the spins of d-quark quasiparticles 
to align the opposite way. The two species being oppositely 
charged, their magnetic moments then add, giving rise to lo- 
cal magnetic fields of order 10'^ G. We estimate the magnetic 
moment density as eH{2S ,i + Sd)/^mc, where S u and Sd are 
the spin densities of u and d quarks respectively. The quark 
mass m is set to g < cr >, g being the Yukawa coupling con- 
stant between the meson and the quark fields, which is set to 
5.4 (see ref. |29]) by fixing the nucleon mass to be 938 MeV. 

The equation of state of this PC matter has been derived and 
presented in ||26[], where we also discuss the structure of neu- 
tron stars that would result when the PC quark matter resides 
in the core. The outer part of the star is made up of nuclear 
matter, for which we adopt the equation of state derived by 
Akmal et al f3^ (APR), used commonly in the literature. The 
boundary between the nuclear and the quark matter PC phase 
is treated as a first order phase transition. The main results are 
as follows. 



A. The quark matter PC core regime 
(750 MeV< <850 MeV) 

1. In this model, the parameter OTq- is constrained to be 
above 750 MeV, to ensure that the two-flavor ground 
state does not fall below the nuclear ground state il8l] . 

2. The maximum mass of neutron stars with PC cores 
works out to be Mmax.pc ~ 1-6 Mq, as in most other 
cases of neutron stars with quark matter core. In con- 
trast, the maximum mass of a purely nuclear star with 
APR equation of state is ~ 2.2 Mq. 

The equation of state for PC quark matter employed by 
us includes gluon interactions as given by Baym 11171 
[stIi . Other prescriptions for this exist in the literature, 
and some attempts have led to the mass limit of a star 
with a quark core being raised to as high as ~ 2 Mq 
MM- 

3. The PC quark matter core occurs only if the mass of the 
star exceeds a threshold value Mj. Below this thresh- 
old mass the star is composed entirely of nuclear matter. 
The threshold mass is a function of the parameter mo-. 
With increasing m^-, the density of phase transition be- 
tween the APR and PC phases rises and so does Mj. 
For an nia- of 815 MeV, Mj works out to be 1.32 Mq. 
A value of nia- > 850 MeV results in Mt rising above 
Mmax.pc and thereby precludes the formation of a PC 
core. Neutron stars with quark matter PC cores can 
therefore exist only if m^r is in the range 750-850 MeV. 

To find the contribution of the PC core to the stellar mag- 
netic field we estimate the magnetic moment fi by integrating 
the magnetic moment density over the volume of the PC core. 
The effective polar field strength at the stellar surface in fully 
relaxed configuration may then be computed as: - 2fi/R^, 
where R is the stellar radius. The result, for two different val- 
ues of mo-, are shown in fig. 1 . As can be seen, has a typical 
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FIG. 1: The polar magnetic field strength at the surface of a neu- 
tron star containing pion-condensed quark matter core. The surface 
field strength is estimated for a fully relaxed dipole configuration 
corresponding to the magnetic moment of the core, shown for two 
values of the parameter nia-, 813 and 815 MeV. The core appears 
above a threshold mass and if present, could generate surface fields 
of strength similar to those seen in magnetars. 

value of ~ 10'^ G, similar to the fields encountered in mag- 
netars. The field strength increases with the mass of the star 
above the threshold mass, and reaches ~ (2 - 3) x lO'^ G near 
the maximum mass. 

If magnetars do indeed owe their magnetic fields to such 
pion cores, their masses must exceed Mj. Ordinary pulsars, 
with much weaker magnetic fields (~ 10'^ G), may then rep- 
resent a neutron star population with masses less than Mj, 
their ma gne tic fields being either inherited from their pro- 
genitors ||33|], or amplified in the crust by thermoelectric pro- 
cesses [34]. This can be tested if the masses of magnetars can 
be measured, perhaps when magnetars in binary systems are 
found. Masses of several pulsars are known, and in most cases 
lie close to ~ 1.4 Mq. 



B. The Purely Nuclear Matter star regime {m^ >850 MeV) 

At high density (above ~ 2-3 times nuclear density Vflie nu- 
clear matter itself can have a neutral pion condensate 0161 [TtIi . 
in which nuclear spins, aligned due to the presence of the neu- 
tral pion condensate, can contribute a strong magnetic field. 
For nicr > 850 MeV the quark matter - nuclear matter transi- 
tion density rises such that only fully nuclear stars are stable 
and can continue to exist. Since it is possible in a purely nu- 
clear matter star at high enough density to have nuclear PC 
core, we then have a nuclear matter analogue for the lower 
threshold mass, M!^ for nuclear magnetars with a new and 
higher maximum mass M^^^^. of the star (as purely nuclear 
EOS's are stifFer). 

We point to this possibility as one object, PSR J075 1 + 1807 
has recently been found to be relatively heavy, with a mass 
of 2.1^Q g Mq |f35l]. While this is still marginally consistent 
with Mmax.pc, if future measurements confirm it to be larger 
than ~ 1.6 Mq then the existence of pure quark matter PC 
core in neutron stars would be ruled out. Core magnetic fields 



through spin alignment may still be produced, but perhaps in 
a pion-condensed nuclear phase |16] or a mixed phase ll36ll 
and the rest of the magnetar scenario would remain the same 
as above. 

This discussion was to indicate that magnetars are possible 
both with a PC quark matter core or with a PC nuclear core. 
However if a neutron star with mass above ~ 1 .6 Mq is found 
then this would preclude pure quark matter core for all stars. 
However, as we have indicated above, it is possible that a dif- 
ferent treatment of gluon interactions may raise the mass limit 
of stars with PC quark cores, as has been obtained in some 
equations of state involving quark matter i3ll l32ll . 



C. Magnetar by birth or accretion 

Neutron stars with a large mass could result either,(i) from 
the core collapse of a rather massive star or (ii) by heavy ac- 
cretion onto a neutron star in a binary system. In either case, 
if the final mass exceeds Mj, a magnetic pion core will form. 
Will one expect to see a magnetar in all such situations? The 
answer depends on the details of the thermal structure in the 
neutron star interior. 

(i) A newly born neutron star in a stellar collapse has a very 
hot interior, facilitating ambipolar diffusion and allowing the 
strong field to emerge at the surface in a short time, as dis- 
cussed above. The more massive the neutron star, the larger 
will be the size of the magnetic core and the quicker will the 
strong field emerge to the surface. It is interesting to note that 
possible magnetar activity has been indirectly inferred in the 
~ 350 yr old supernova remnant Casseiopeia A [Issll, which 
originated in the explosion of a star as massive as ~ 25 Mq 
[39]. Most importantly, recent observations |40] have uncov- 
ered several cases of association between magnetars and mas- 
sive star clusters indicating that the progenitors of these neu- 
tron stars were more massive than ~ 30 Mq. A similar con- 
clusion is drawn in the case of AXP IE 1048. 1-5937 which is 
associated with a stellar wind-blown shell of neutral hydrogen 

(ii) A neutron star accumulating matter via accretion, on 
the other hand, is old and its interior is relatively cold. Heat- 
ing due to the accretion process is not expected to raise the 
interior temperature above ~ 10^ K |41]. The extreme tem- 
perature sensitivity of the ambipolar diffusion rate will then 
delay the emergence of the field at the surface, perhaps for 
such a long time that the magnetar property would never be 
visible. This may be the reason why the surface magnetic field 
PSR J075 1-1-1 807 remains low II35I1 despite its mass growing 
to a large value. 



IV. CONCLUSIONS 

We have found that in certain allowed equations of state 
of neutron stars a strongly magnetized pion condensed core 
occurs at high density, for neutron star masses near the maxi- 
mum mass. We have argued that such stars would have prop- 
erties remarkably similar to the observed magnetars. The 
magnetic field of the core would initially be squeezed deep 
into the star due to screening currents set up in the material 
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surrounding the core. The free energy in this field configura- 
tion would be released over a time scale of ~ 10"^ yr and power 
the magnetar activity. Along with this, the magnetic field at 
the surface will rise, possibly turning Anomalous X-ray Pul- 
sars into Soft Gamma Repeaters. If this model for magnetars 
is correct, then magnetars would be expected to have masses 
higher than a certain threshold, and would therefore be ex- 
pected to be associated with supernova explosions from rela- 
tively massive stars. Neutron stars with magnetic cores may 
also be produced by heavy accretion onto normal neutron stars 
in binary systems, although it seems less likely that the mag- 
netic field will emerge to the surface quickly enough to exhibit 
magnetar activity in this case. If a magnetar in a binary system 
is found, several of these predictions could be tested. 
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